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Abstract: We give necessary and sufficient conditions for two sub-vectors of a random 
vector witli a multivariate extreme value distribution, corresponding to the limit distribu- 
tion of the maximum of a multidimensional stationary sequence with extremal index, to 
be independent or totally dependent. Those conditions involve first relations between the 
multivariate extremal indexes of the sequences and secondly a coefficient that measure the 
strength of dependence between both sub-vectors. The main results are illustrated with 
an auto-regressive sequence and a 3-dependent sequence. 

1 Introduction 

Multivariate Extreme Value Analysis is frequently applied in the context of modeling 
environmental data, for which the phenomenon of dependence is often intrinsic. This 
paper focuses on the characterization of total dependence and of independence of two 
multivariate extreme value distributions. 

Let X = {X^'^'* = {Xn,i, ■ ■ ■ , Xn^d}n>i be a d— dimensional stationary sequence with 
common distribution function (d.f.) Q{x^'^^) = Q{xi, . . . ,Xd), x^'^) E M*^, and M„ = 
{Mn^i, . . . , Mn^d) the vector of pointwise maxima, where M„^j is the maximum of i-th 
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component of X. Denote M„ = (M„_i, . . . ,Mn,d) the corresponding vector of pointwise 
maxima of the associated d-dimensional sequence, X = {X^'^^}„>i, of independent and 
identically distributed (i.i.d.) random vectors having the same distribution function Q. 

In this multivariate setting operations among vectors are defined componentwise, that 
is for each d > 1 and a^*^) , b^*^) G W^, a^'^) < b^*^) , if and only if aj < bj , for all j = 1, 2, . . . , d. 

If there exist sequences {{an,i > 0, . . . , an,d > 0)}„>i and {(6n,i, • • • , bn,d)}n>i, such 
that for u(xW) = {u„(x('^)) = {an,ixi + bn,i, an,dXd + bn,d)}n>l, 

P (m„ < u„(xW)) = {^".J- ^ + ^".i} j ^(^^''^)' ^^"^ ^ 

where G is a d.f. with non-degenerate margins, then Q is said to be in the max-domain 
of attraction of G (Q G D{G)) and G is said to be a Multivariate Extreme Value (MEV) 
distribution function. 

We will assume, without loss of generality, that the univariate marginal distributions 
of G are equal to F. 

It is well known that the relationship between the d.f. G(x('^)), x^*^) G W^, and its 
marginal distributions F{xj),j = can be characterized by its copula function, 

Dq, which exhibits a number of interesting properties, namely its stability equation: 

Dhiyi. • • • , = Dcivl . . . , y'd), > and (yi, . . . , G [0, 1]"' . (1.1) 

If the stationary sequence X satisfies some mixing conditions, £)(u„(x('^))) of Hsing 
([4]) or A(u„(xW)) of Nandagopalan ([6]), and 

where is a d.f. with non-degenerate components, then R is also a MEV d.f.. The 
MEV d.f a and G can be related through the multivariate extremal index function, 
6{t^'^) = 6 {ti, . . . ,Td) introduced by Nandagopalan ([6]), which is a measure of clustering 
among the extreme values of a multivariate stationary sequence. 

Definition 1.1 A d— dimensional stationary sequence X, is said to have multivariate 
extremal index 0-^{t'''^^) G [0, 1] , if for each t^'^^ = (ti, ... ,Td) G M^|_ there exists u^^^ = 
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(Vi ' • • ■ ' )' ^ - ^' satisfying 

nP{Xij > u^^^h > Tj, j = l,2,...,d, P (m„ < n^;:'"')) > G{t^<''^) and 

As in one dimension, the extremal index is a key parameter when relating the properties 
of extreme values of a stationary sequence to those of independent random vectors from 
the same c/— dimensional marginal distribution. However, unlike the one dimensional case, 
it is not a constant for the whole process, but instead depends on the vector t^'^\ 

It is now clear that the existence of 6^(t^'^^) allows us to write 

ii'(x('^)) = G^''(^''')(x('^)) with Tj = Tj{Xj) = -\ogF{xj), j = l,2,...,d. 

Taking d = p + q, it follows, as a consequence of the definition of multivariate extremal 
index, that the sequences X(?') = {X^f^ = {Xn,i, . . . , Xn,p)}n>i and X(«) = {X^*^ = 
{Xn^p+i, . . . , Xn,p+q)}n>i have, respectively, extremal indexes 

^x(rt^^(p)^^ lim e^{T^+i^) and ^^^'^(t^^)) = lim e^{T^+i^). 

-rj-^0+ Tj->0+ 
j=p+l,...,p+q j = l,...,p 

In the notation of the extremal index we shall omit the sequence, whenever it is clear 
by the context and the argument of the function. 

Hereinafter, let Y = {Yi, . . . ,Yp+q) and Y = {Yi, . . . ,Yp+q) be, respectively, 
two random vectors with distribution functions Q^^i'^'-^'^"^) and G, where Y^^*) = (Yi, . . . , Yp) 
and Y^^) = (l^+i, . . . , denote two sub- vectors of Y and Y^*') = (?!,..., Yp) and 

Y^*) = (1^+1, . . . , Yp+q) two sub- vectors of Y. 

In section 2 we discuss conditions under which Y*^) and Y^^) are independent or 
totally dependent. These conditions are established first by relations between the extremal 
indexes ^^(t(p+9)), e^''''(T(f)) and ^^'''(t^'?)) and secondly by a coefficient that measures 
the strength of dependence between Y*^) and Y^^). 

The main results are illustrated in section 3 with an auto-regressive sequence and a 
3-dependent sequence. 
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2 Main results 



If the d.f. Q belongs to the domain of attraction of a MEV distribution, G, and X has 
extremal index 6I(t(p+«)), t(p+«) = ( Ti, . . . , Tp+q) G , then we have 

Ge(r(''))(xW)G^(-(''^)(xW) < G^(-^''+^')(x(f+«)) < min{G^(-'^')(x(P)),G^(-'^')(x(«))}, 

(2.2) 

for each x(p+«) G W+i and r,- = - logF(xj), j = 1, . . . ,p + g. 

The inequality on the right holds true for every multivariate distribution, while the 
inequality on the left is a property of MEV distributions. The lower bound corresponds 
to the case where Y^p^ and Y(«) are independent and the upper bound corresponds to the 
case where Y^p) and Y(«) are totally dependent. 

From (2.2) we obtain the following bounds for the multivariate extremal index function 

max{g(TW)7(TW),g(T('^))7(TW)} g(r(p))7(rW) + g(T(^))7(Tfa)) 

7(T(f+9)) - ^ ^ - 7(T(f+9)) ' ^ ' 

where 

^(r(P+«)) = -logG(F-i(e-"i), . . . , F-i(e-"f+^)) = lim nP fx5^+^) ^ 4^^'^'^^) > 
7(t(p)) = lim 7(T(f+«)) and 7(r(«)) = lim 7(t(*'+«)). 

Tj-i-0+ 7-^-^0 + 

i=P+i,...,P+c3 j=i,...,p 

The next result follows naturally from these bounds. 

Proposition 2.1 Suppose that Q G D{G) and X /tas extremal index ^(r^+^)), t^+^) G 

(i) /f and Y^^^ are independent, then Y^^ and Y^^^ are independent if and only if 

^(^^^^)7(r(^)) + g(T(-'))7(Tfa)) 
^" 7(T(f))+7(TW) ' + ■ 

(ii) // Y^P) and Y^'^) are totally dependent, then Y^^ and Y^''^ are totally dependent if 
and only if 

_ max{g(rW)7(T(^)),g(r('^))7(r(^))} ^ 
^ ^ max{7(T(f)),7(T('^))} ' 
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The necessary and sufficient conditions for Y^^^ and Y^^^ to be independent or totally 
dependent given in the previous result demand the evaluation of the extremal index func- 
tion 0(t(^+^)), in each point t^p+i") g M^"*"^. Nevertheless this task can be simplified with 
the characterizations, given in Ferreira ([3]), for independence and total dependence of 
the multivariate marginals of a MEV distribution. These characterizations are essential 
to prove the following propositions which guarantee that the independence or total de- 
pendence between Y^*') and Y^^) only depends on the value of the extremal index in some 
points. 

Proposition 2.2 Suppose that Q G D{G) and the sequence X = {X^^'^^}„>i has extre- 
mal index 0(t(p+«)), t(p+«) € 

The sub-vectors and Y^'^) are independent if and only if 



' ' 7(l(p+9)) 
where 1^*^^ = (1, . . . , 1), A; > 1, denotes the k-dimensional unitary vector. 



(2.4) 



Proof: Suppose that Y^^ and Y^^^ are independent. Since (2.3) holds for all 
riP+i) e MF^'', we have in particular for t^+i^ = {t,...,t) G M^+«, with 
T = t{x) = - \ogF{x), X eR, 

^ 0(rW)7(TW) + g(Tto))7(Tfa)) 

Now attending to the fact that 9{ct^'^^) = 9{t^^^) for each r^''^ € M^, /c > 1 and c > 0, we 
can write 

e{T^P+i^) = e(l(P+«)), 0(t(p)) = ^(1^^)), ^(r^")) = ^(i(«)), 
and from (1.1), for all t(p+«) = (r, . . . ,t) G M^+^ 

7(r(f+'?)) = -\ogG{F-\e-^),...,F-\e-^)) = -\ogDG{e-\...,e-^) 

= -logD£(e-i,...,e-i) = T7(l(^+''^ (2.5) 

j{t^p)) = Tj{l^P^) and 7(t(«)) = r7(l(«)). Equality (2.4) is now straightforward. 
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On the other hand if (2.4) is verified, then for x^^^^^ = [x,. . . ,x) we have 

= exp(-r(0(l(f))7(l^)) + ^(l(^^)7(l(^^))) = GY(.)(x(^))GY(.)(x(«)) 
and from Proposition 2.1 ([3]) we conclude that Y^^^^ and Y^^) are independent. □ 
Proposition 2.3 Suppose that Q G D{G), X = {X^f+^^} n>i has extremal index 

(i) // Y^*") and Y^^) are totally dependent then there exists t'^+i^ g M^"*"* with 
Tj = Tj(xj) = — logF{xj), Xj G M, j = 1, . . . ,p + q, such that 

7(T(f))0(T(f)) = 7(T(«))e(T(«)) = 9iTi... = 9p+,Tp+g = d>0 
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and 0(t(p+'?)) = 7 

(ii) // there exists t(p+'?) G M^J."^^ mi/i = Tj(a;j) = - log F{xj), xj eR, j = 1, ...,p + 
such that 

7(T(f+^))^(T(f+^)) = ^in . . . = ep+5Tp+5 = d > 0, 
then Y(p) anc? Y^*) are totally dependent. 

Proof: (i) From Proposition 2.1 ([3]), if Y^^) and Y^''^ are totally dependent, then there 
exists T^+'^^ G such that 

9{T^^h{T^^) = 0(t(«))7(t(«)) = 0(T(^'+«))7(T(f+«)) = d = 9in = ...= 9p+,Tp+„ 

with d g]0, 1[. Hence 



-log£)G(exp(-Ti), . . . ,exp(-rp+g)) 
1 1 



■logZ^G (cxp(-^),...,exp (-li^)) 
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-{p+'i) 



□ 
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Another way to look at issues concerning independence or total dependence is to 
use parameters that measure the strength of dependence between Y^*') and Y^^). We 
therefore define, in the following result, the dependence structure of Y^*') and Y^^) through 
the coefficient e{Y(^\Y(«)) q£ pgrreira ([3]). This coefficient emerged from the extremal 
coefHcient of Y, e^, defined in Martins and Ferreira ([5]) as 

and the relationship 

It is then defined as £(y(p',y(9)) _ ^^(j i^^g following interesting properties. 

Proposition 2.4 

^ ' ~ (?(1(P))7(1(P))+0(1('J))7(1(9)) 

(ii) g(Y<''',Y('2)) _ ^ if and only ifY^P^ anclY^''^ are independent. 

(iii) // Y^*') and Y^'^) are totally dependent, then g(Y'''',Y('j)) _ — ,e } 

Proof: (i) Since 

^(yM,yM) ^ g(i(p+^))iogG(x(P+g)) 

log G(x(p)) + 6'(l(«)) log G(x(3))) 

the result follows from (2.5). 

(ii) It is an immediate consequence of (i) and Proposition 2. 

(iii) If Y^^) and Y^'^^ are totally dependent then from (2.2), we have 

(y(p),yW) _ -logmmlG"'-' '(x''"), G"'-' '(x^^'^)} _ maxjf } 
^ ' ~ ^(iW) log G(x(p)) + 0(1(9)) log G(x(9)) ~ eY(p)+gY(«) 

□ 
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3 Examples 



Example 3.1 Let {i^}n>i be a sequence of i.i.d. random variables with common d.f. F 
and consider an auto-regressive sequence of maxima {X„}„>i defined by 

Xn = max{y„,y„+i}, n > 1, 

with marginal distribution function F^. 

Let {un''^}n>i, i = 1, . . . ,p, and {vn^^}n>i,j = p + I, . . . ,p + q, be sequences of real 
numbers such that lim„^oo^(l ~ F'^{un'^)) = ti and lim^^oo ^^i^^)) = t'-. 

The sequences {Xn)n>i and {—Xn}n>i have, respectively, extremal indexes 9\ = 1/2 
and 62 = 1. 

T-> -viP+l) J ^n,i = Xn , i = l,...p (p) 

For sequences X^^^^'' = < , X^^^ = 

[Xn,i = -Xn , i=p+l,...,p + q 

and Xn*-* = {—Xn, ■ ■ ■ , —Xn), we have 

lim P(M(f) < 4-^''')) = expf-i m^r.y 

hm P(M(f) < ui-'"')) = exp f- mc^r,) , 

hm P fM^"^) < . . . , -yp+^M = exp ( - max r' ) . 

Since the order statistics maximum and minimum are asymptotically independent ([2], [7]) 
we obtain 

P{Mn < (4^1), . . . , u'^:''\v'a^+'\ > exp (-1 max r, - max rj) 

and consequently j{t^~^'^^) = j{t^^) + 7(t('^) = maxr^ + max t', and 
6l(r(P+9))7(T(P+«)) = i maxr,- + max r'. Therefore 

l<j<p P+l<j<P+9 



7(t(p)) + 7(t(9)) 

Example 3.2 Let U = {Un}n>i be a sequence of i.i.d. random variables with common 
d.f. H in the domain of attraction of the extreme value distribution F, and independent 
of the i.i.d. chain J = {Jn}n>i such that P(Ji = 0) = P(Ji = 1) = 1/2. 
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Let us consider a stationary 1-dependent sequence Z = {Zn}n>i, defined as Zn = Un 
if J„ = and Zn = Un+i otherwise, and let v = {vn}n>i be a sequence of normalized 
levels to Z, and consequently also to U. 

We can now define a 3-dependent stationary sequence X = {X„ = {Xn,i, Xn,2, Xn,3)} 

as 

{Xn,l,Xn,2, Xn,3) = {Zn, Zn+2, Zn+l), n > 1, 
with common distribution function 

1 1 1 

T{xi,X2,X3) = ^YiH{Xi) + -H{xi)H{mm{x2,X3}) + -H{x2)H{mm{xi,X3}) 

i=l 

belonging to the domain of attraction of 



G{X1,X2,X3) = < 



F{xi)F{x2)F2{x3) , Xi<X3AX2<X3 

F(xi)F3(x2)-F4 (xa) , Xi<X3AX3<X2 

F2(xi)F(x2)-F^(x3) , X3<XiAX2<X3 

3 3 

F4(xi)F4(x2)F(x3) , X3<XiAX3<X2 



Now applying Proposition 2.1 of Smith and Weissman (1996) ([8]) to the sequence 
= {maix{Xn,i,Xn,2,Xn,3} = max{Z„, .Z„+2}}n>i which verifies the condition 
D^^\vn), A; = 2, of Chernick et al. ([1]), we easily obtain 

e^il 1 I) = lim ^ ^""^ "^^^ > > ma^ {Z2, Z3, Z4}) ^ _3_ 

^' ' ^ P{iaax{Zi,Z2,Z3} > Vn) 10' 



For random vectors Y = {Yi,Y2,Y3j and Y = {Yi,Y2,Y3) with d.f. = G and 
Gy = Y(2) = {Y^,Y2) and Y(i) = Y3 we obtain = |, = e^C^) ^ 2, 

= 1, = 3 ^ ^Y(^). Consequently = 1, = | and from 

Proposition 4 we conclude that neither Y^^) and Y^^^ nor Y^^) and Y^^^ are independent. 

Nevertheless, there exists r^^^ = (1,1,1) such that ^{t^-^^)0{t^'^^) = OiTi = 62T2 = 
^3''"3 = f s-iid attending to Proposition 3 we can say that Y^^^ and Y^^^ are totally 
dependent. 
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